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Introduction and Background 


24.08.2009 
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COMMON WISDOM ON LSB 


Ml 


./A] jy 


* Separation under adverse pressure gradient 

* Reversed flow region enclosed ^ Shear layer 

* Kelvin-Helmholtz instability ^ transition to turbulence 

* Turbulent strong mixing forces re-attachment 
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(Horton 1968) 
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• Two different classes of eigenmodes: 

Shear-layer/KH Global mode 



24.08.2009 Th, Hein & Dallmann Phil. Trans. R. Soc. Lond. A (2000) 358: 3229-3246 
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LPT FLOWS 

- THE (SL) WAKE MODE 



• Amplification of the wake mode gives rise to the 2D time- 
periodic basic state 



qO, y, t) = q(>, y) + £ qO, y ) exp © 


24.08.2009 


Abdessemed, Sherwin, Th, JFM (2009) 628: 57-83 
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LPT FLOWS - THE BUBBLE MODE 



■ 




The bubble mode is 

stronger damped than the wake mode is 


disturbance vorticity spanwise velocity 
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Abdessemed, Sherwin, Th, JFM (2009) 628: 57-83 
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A UNIFYING PERSPECTIVE 


THE GLOBAL MODE OF LSB 



• The struture of the amplitude functions 

of the minimally damped global eigenmode of the LPT blade... 

spanwise w(x,y) 



...is analogous to that on the flat-plate and the NACA 0012 aerofoil 


24.08.2009 
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GLOBAL INSTABILITY OF LSB 


IS 


Configuration 

Authors 

Methodology 

Airfoils 

Gaster (1964) 

- 

APG flat plate BL 

Hammond & Redekopp (1998) 

Absolute/Convective 
instability analysis 

-II- 

Theofilis (1998, 1999) 
Theofilis, Hein & Dallmann (2000) 

BiGlobal 

instability analysis 

Backward- and forward-facing 
steps 

Barkley, Gomes, Henderson (2001) 
Marquet et al. (2006) 

Marino & Luchini (2009) 

BiGlobal 

instability analysis, 
DNS 

Low Pressure Turbine blades 

Abdessemed, Sherwin & Theofilis 
(2004; 2006; 2009) 

BiGlobal analysis, DNS 

Shock / Boundary-Layer 
Interaction 

Boin, Robinet (2004; 2007) 

BiGlobal instability analysis, DNS 

NACA Airfoils at AoA 

Theofilis, Barkley & Sherwin (2002) 
Crouch et al. (2007; 2009) 

BiGlobal instability analysis, 
DNS 

Rounded open cavity, 
S-shaped duct 

Hoepffner et al. (2007), Akervik et al. 
(2007), Marquet et al. (2008) 

BiGlobal instability analysis, DNS 

Bump-induced separation 

Ehrenstein et al. (2007) 

BiGlobal instability analysis, DNS 

Low Pressure Turbine Cascade 

Brear & Hodson (2004) 

Experiment 

Iced Airfoil 

Jacobs & Bragg (2006) 

Experiment 
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GLOBAL INSTABILITY & FLOW TOPOLOGY - 1998 


• Dallmann’s conjecture on the topological flowfield changes due 
to the global mode 



24.08.2009 


Th, Hein & Dallmann Phil. Trans. R. Soc. Lond. A (2000) 358: 3229-3246 
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GLOBAL INSTABILITY & FLOW TOPOLOGY - 1999 






• Amplification of the global mode leads to three- 

dimensionalization of the primary reattachment line in 3d... 



...while it leaves the primary separation line unaffected 


24.08.2009 


Th, IUTAM Laminar-turbulent transition V, Sedona (1999) 663-668 
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The present contributions 


24.08.2009 
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OUTLINE 


Incompressible Flow 

• Some numerics on the direct/adjoint EVP 

• U-separation on a flat-plate 

• Stall cells on a NACA0015 airfoil 

• Receptivity and Sensitivity of LSB on a flat plate 
Compressible Flow 

• The Howarth/Briley bubble in subsonic flow 

• Shock/Boundary-Layer Interaction (SBLI) on a semi-infinite plate 
at Mach 2 

• SBLI on a 45° wedge 

Focus on the (stationary, self-excited) global mode alone 


24.08.2009 
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Numerical considerations 


24.08.2009 
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MASSIVELY PARALLEL SOLUTION OF THE EVP 


• Ten years later. 


Incoming TS waves at the inflow boundary 

(besides the homogeneous Dirichlet BCs used in THD) 

- 400x100 points per spatial direction used to resolve the 
amplitude function of the eigenmodes 

(as opposed to 60x40 in THD) 

- 25k CPU hrs ( = 25 hrs x 1000 procs) per wavenumber 

(as opposed to 24 CPU hrs serial computing time on a single 
processor in THD) 

...permit solving substantially larger (but not all) problems, 

• on Mare Nostrum at BSC and (the Petaflop) Blue-Gene/P at FZJ 

(but Jurgen Seidel finally got access to DoD HPC machines) 


24.08.2009 


Rodriguez & Th AIAA J (2009 - accepted) 
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Classification of global modes of LSB flow 


24.08.2009 
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GLOBAL MODES IN LSB FLOW ON A FLAT PLATE 


Case: Re 8 ,=500 at inflow, Re g ,=700 at outflow. Separation Bubble 
Boundary conditions: Dirichlet at inflow & Extrapolation at outflow 


- 0.02 
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Walll Modes 
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24.08.2009 


Rodriguez &Th AIAA 2008-4148 
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TYPICAL FLAT-PLATE GLOBAL EIGENSPECTRUM 


Composed of discrete and continuous branches 

Focus discussion on two traveling ( and th e l eading stationary) 
global eigenmodes 

i i i i i i i 



A 

* 





/•"A 

0.02 0.04 0.06 0^8 0.1 

Rg (a>) 


0.12 0.14 0.16 


24.08.2009 


Rodrfguez & Th AIAA Pap 2008-4148 
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Comparison of traveling modes with OSE 


• Extract “spatial amplification” from successive peaks of the (2d) 
global eigenfunction; 

• Run a spatial OSE at successive downstream locations; 

• Compare 


10 ' 

£ 

< 


10 ' 
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Extracting physics out of eigenspectra: 
• U -separation on a flat plate 


24.08.2009 
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U-SEPARATION 


Fluid Mechanics Group -- MfrME - Windows Internet Explorer 



T H_, http://www.mame. mu, oz.au/fluids/ 

V \*t X 

Google | P | T | 

^Convert t ^Select 



^ $$ Fluid Mechanics Group -- M&ME 

a - m 

E^Q T [T^ Pagina ^ Herramientas ^ 


Fluid Mechanics Group 

Walter Bassett Aerodynamics Laboratory 
Department of Mechanical Engineering;, University of Melbourne, Australia 



The above digram represents a local Taylor series expansion solution oftheXavier-Stokes equation for a simple U-separation. 
(F rom AT. Pern - & MS. Chong. Annual Review of Fluid Mechanics. 1987. YoL 19. pp.125-55). 

Academic Staff 
Research Staff 


0 
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“U-SHAPED SEPARATION” - 1984 



FJG, J 7 : of the fret si ream surface sheets and 

WLilJ streamlines of si nn pile U-sliapcd separation 


24.08.2009 


Hornung, Perry Z. Flugwiss. Weltraumforsch.(1984) 8(2):77-87 
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For example, Fig. 6 shows a three-dimensional separation 
pattern that has been classified as a U separation ( Perry and 
Hornung 1 T ). The flow field has been obtained by solving the 
Navier-Stokes equations locally (see Perry and Chong 18 ) 
using a third-order Taylor series expansion, and assuming 
that the flow pattern is symmetrical (ie,, x, and ± 3 are as- 


24.08.2009 


Chong, Perry, Cantwell Phys. Fluids A (1990) 2(5):765-777 
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TYPICAL FLAT-PLATE GLOBAL EIGENSPECTRUM 


Composed of discrete and continuous branches 

Focus on ( two traveling and ) the leading stationary global 
eigenmode 

i i i i i i i 



A 

* 





/•"A 

0.02 0.04 0.06 0^8 0.1 

Rg (a>) 


0.12 0.14 0.16 


24.08.2009 


Rodrfguez & Th AIAA Pap 2008-4148 
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Critical Points of ( BF + e A ) 


Reconstruct a composite flowfield, using the steady basic state 
and the leading global eigenmode 

Identify critical points and their nature 


220 


260 220 240 260 220 240 260 220 240 260 

XXX 


200 


250 


300 


350 


S: Saddle N: Node F: Focus C: Center 
(s): stable (u): unstable 


24 . 08.2009 
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Critical Points of ( BF + e S ) 





• Critical point analysis of composite flow 

Eigenvalue spectrum BF + e (Global Mode) 



24.08.2009 


Rodriguez & Th J Fluid Mech (2009 - submitted) 
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Present computations 



Bippes (198?) 


24.08.2009 
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Extracting physics out of eigenspectra: 

• Stall Cells on an Airfoil 


24.08.2009 
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24.08.2009 


Bippes (198?) 


- courtesy of H. Hornung via H. Fasel 


STALL CELLS ON AIRFOILS 




NACA 4415 


A = 3,1 
a = 21,5° 
Re= 2,1x10 6 



Wandstromlinien auf der Saugseite des Tragfliigels 

DFVLR 

v pc 

bei teilweise abgeloster Stromung 
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STALLED AIRFOIL 


Critical point analysis of composite flow 


Eigenvalue spectrum 


Leading stationary eigenmode 



u-perturbation 


pressure perturbation 



24.08.2009 


Kitsios, Rodriguez, Th, Ooi, Soria J Comput Phys 228:7181-7196 (2009) 
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Critical Points of ( BF + e GlobalMode ) 



• Stall-cells appear on the airfoil, as result of linear modal 
amplification of the leading stationary global mode 



24.08.2009 


Rodriguez & Th Theor Comp Fluid Dyn (2009 - in review) 
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Flow Control of LSB 


24.08.2009 
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• The kernel of the direct/adjoint coupling (Luchini 2003) 
permits identifying the region of maximal feedback 
region... 



• ... it coincides with the primary recirculation center 

• (a result not unexpected, really) 


24.08.2009 


Rodriguez, Valero, Ekaterinaris, Th IUTAM Laminar-Turbulent Transition 

Stockholm (2009) 
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RECEPTIVITY OF LSB ON A FLAT-PLATE 


• Solution of the adjoint eigenvalue problem permits identifying 
the maximally receptive region (dominant streamwise 
velocity perturbation) 



• ... it coincides with the primary separation line 

• NB. Analysis is meaningful due to homogeneous Dirichlet BCs 
imposed upstream (direct) / downstream (adjoint) EVP 


24.08.2009 


Rodriguez, Valero, Ekaterinaris, Th IUTAM Laminar-Turbulent Transition 

Stockholm (2009) 


33 




NASA/CP — 20 10-216112 238 



SCHOOL OF AERONAUTICS 


UNIVERSIDAD POLITECNfCA DE MADRID 


WHAT CAN THIS KNOWLEDGE DO FOR US? 


• Predict the actuator placement without prior knowledge (or 
ad-hoc assumptions ) 



0 . 


I 

/ 

/ 
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Hill AIAA 1992-0067; Giannetti, Luchini JFM (2007) 
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Compressibility effects 


24.08.2009 


UNIVERSIDAD POLITECNfCA DE MADRID 


- subsonic flow 
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MACH # EFFECT ON HOWARTH/BRILEY LSB 


• Linear pressure increase (compressible) v linear deceleration 
(incompressible) in the tree-stream 
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• No major differences in integral quantities have been found, ie. 

• Quantitative but no qualitative differences of instability results 

• Work in progress, in collaboration with Stanford U. 


24.08.2009 


Valero & Th AIAA Pap 2008-0593 
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Compressibility effects - supersonic flow 


24.08.2009 
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SHOCK / BOUNDARY-LAYER INTERACTION 


• Typical LES result at Re s = 10k 



• Origin of unsteadiness? 

• Origin of spanwise periodic structures ? 


24 . 08.2009 


Ekaterinaris AIAA J (2004) 42(3):493-500 
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SBLI & 45° COMPRESSION CORNER - BFs 






• Basic states obtained by high-order accurate, low-diffusion 
WENO schemes (Ekaterinaris) 

• Flat plate leading edge included at upstream end of both 
domains 

• Finite-angle wedge in compression corner (no triple-deck) 

• Both flows at 5000 < Re L < 20000, Ma = 2 and 2.5 

• Unsteadiness identified by DNS at Re L » 25k 



24.08.2009 


Ekaterinaris AIAA Pap. 2008-4172 
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• Compressible BiGlobal EVP solved at 5000 < Re L < 20000 by 

• Spectral collocation (SBLI on flat plate) 

• Immersed boundary (45° compression corner) 

• Entire BF domain considered - shocks not treated ( t Crouch ) 

• Qualitatively reasonable results,... 



• ...but spectra are not converged (despite very large resolutions) 

• Work in progress, in collaboration with with Patras U. 


24.08.2009 


Rodriguez, Valero, Ekaterinaris, Th IUTAM Laminar-Turbulent Transition 

Stockholm (2009) 
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SUMMARY ( I / II ) 


• In the last 3 years, global linear instability of LSB has been 
revisited, using state-of-the-art hardware and algorithms 

• Eigenspectra of LSB flows have been understood and classified 
in branches of known and newly-discovered eigenmodes 

• Major achievements: 

• World-largest numerical solutions of global eigenvalue 
problems are routinely performed 

• Key aerodynamic phenomena have been explained via 
critical point theory, applied to our global mode results 

• Theoretical foundation for control of LSB flows has been laid 


24.08.2009 
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• Global mode of LSB at the origin of observable phenomena 

• U-separation on semi-infinite plate 

• Stall cells on (stalled) airfoil 

• Receptivity /Sensitivity /AFC feasible (practical?) via 

• Adjoint EVP solution 

• Direct/adjoint coupling (the Crete connection) 

• Minor effect of compressibility on global instability in the 
subsonic compressible regime 

• Global instability analysis of LSB in realistic supersonic flows 
apparently quite some way down the horizon 

(Saric joke # 431 : “an imaginary, ever-receding line”) 


24.08.2009 
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FUTURE CHALLENGES AND WORK DIRECTIONS 


• From a physical (and theoretical point of view 

• Continue with the analysis of EVP results 

• Instability analysis of turbulent flows (obtain EVP spectra of 
2.5D flows) 

• Embed instability analysis into (AFC) flow control concepts 
(for 2.5D flows) 

• Extend global instability tools in 3D flows 

• From a numerical point of view 

• Devise new algorithms for the efficient parallel solution of 
(massive) eigenvalue problems 


24.08.2009 
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Thank You 


24.08.2009 
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GLOBAL FLOW INSTABILITY AND CONTROL IV 



Crete, Greece, Sept 28 - Oct 2, 2009 
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ADDITIONAL MATERIAL 
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BASIC LSB FLOW ON A FLAT PLATE 



Boundary Layer transformation: 


5 L 


n = y 


u. 


VX 


^ = y/U^vxf(Z,Tj) 


Separated states recovery: 

- Reyhner and Flugge-Lotz approximation 

- Displacement thickness imposed 
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COMPARISON WITH TEMPORAL OSE 


Basic Flow: Artificial parallel Blasius 

Mack’s Case: Re = 580, a = 0.1 79 (Mack jfm 1976) 

Streamwise extension: L =10 x 2n / a 
Periodic boundary conditions 
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FURTHER COMPARISONS WITH TEMPORAL OSE 


Basic Flow: Artificial parallel Blasius 

Mack’s Case: Re = 580, a = 0.1 79 (Mack jfm 1976) 

Streamwise extension: L =10 x 2n / a 


Boundary conditions: 

X Periodicity - Fourier 

+ Periodicity - 
Chebyshev 

O Robin inflow & outflow 

n Robin inflow & 
Extrapolation outflow 

''‘"- 0.1 - 0.08 - 0.06 - 0.04 - 0.02 0 0.02 0.04 0.06 0.08 0.1 
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Case: Re s ,=500 at inflow, Re 5 .=700 at outflow 

Analysis 1 : 

- Basic Flow: Artificial parallel Blasius 

- Boundary conditions: Robin* at inflow & outflow 

Analysis 2: 

- Basic Flow: Artificial parallel Blasius 

- Boundary conditions: Robin* at inflow & extrapolation at outflow 

Analysis 3: 

- Basic Flow: Real Blasius boundary layer 

- Boundary conditions: Robin* at inflow & extrapolation at outflow 

* Robin boundary condition of Ehrenstein & Gallaire JFM 2005: 
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FURTHER COMPARISONS WITH SPATIAL OSE 



Case: Re 5 *=500 at inflow, Re 5 *=700 at outflow 

Robin boundary condition evaluated at: oo 0 = 0.13 


Analysis: 


1 X Parallel Blasius, 
Robin + Robin 


2 + Parallel Blasius, 
Robin + Extrapolation 


3 O Real Blasius, 
Robin + Extrapolation 

0 0.05 0.1 0.15 0.2 0.25 0.3 


0.02 1— r 



- 0.02 


-0.04 


-0.06 


-0.08 


- 0.12 


-0.14 
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